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Abstract The recently developed void expansion 
method (VEM) allows for an efficient generation of 
porous packings of spherical particles over a wide range 
of volume fractions. The method is based on a ran- 
dom placement of the structural particles under addi- 
tion of much smaller "void-particles" whose radii are re- 
peatedly increased during the void expansion. Thereby, 
they rearrange the structural particles until formation 
of a dense particle packing and introduce local hetero- 
geneities in the structure. In this paper, microstructures 
with volume fractions between 0.4 and 0.6 produced by 
VEM are analyzed with respect to their degree of het- 
erogeneity (DOH). In particular, the influence of the 
void- to structural particle number ratio, which consti- 
tutes a principal VEM-parameter, on the DOH is stud- 
ied. The DOH is quantified using the pore size distribu- 
tion, the Voronoi volume distribution and the density- 
fluctuation method in conjunction with fit functions or 
integral measures. This analysis has revealed that for 
volume fractions between 0.4 and 0.55 the void-particle 
number allows for a quasi-continuous adjustment of the 
DOH. Additionally, the DOH-range of VEM-generated 
microstructures with a volume fraction of 0.4 is com- 
pared to the range covered by microstructures gener- 
ated using previous Brownian dynamics simulations, 
which represent the structure of coagulated colloidal 
suspensions. Both sets of microstructures cover simi- 
larly broad and overlapping DOH-ranges, which allows 
concluding that VEM is an efficient method to stochas- 
tically reproduce colloidal microstructures with varying 
DOH. 
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1 Introduction 

The mechanical properties of colloidal particle gels are 
of great importance in nature and in many of today's 
technologies. Sediments pQ, clay [2], some food [3] or 
paints and coatings [4] are but a few examples. Typ- 
ically, colloidal gels exhibit very complex mechanical 
behaviors such as shear thickening [5] , thixotropy [5] or 
aging 7 . From a scientific perspective, these mechan- 
ical properties are often studied using random sphere 
packings. They provide ubiquitous model systems for 
colloids - or granular materials in general - and allow 
for a systematic study of the mechanical properties as a 
function of the various parameters as, for example, the 
volume fraction [8], the particle size distribution [9], 
material properties such as the particles' friction coef- 
ficients [TU] or adhesive forces [TTj . 

Our present research focuses on the influence of the 
microstructural arrangement of the primary particles 
on the structure's mechanical properties. This influence 
is often observed implicitly in experimental and simu- 
lated mechanical tests on gels or particle packings dif- 
fering in preparation history. Macroscopic stress pro- 
files, for example, were found to depend strongly on 
the sample preparation procedure and thus on the mi- 
crostructure [T2]. Franks et at [13] have investigated 
the mechanical properties of sedimented aggregates of 
submicron alumina particles and, in particular, the in- 
fluence of the aggregate size (i.e. the microstructure) on 
the rheological properties. They found that sediments 
formed from larger aggregates exhibit higher shear and 
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compressive yield strengths. In [HlfTMTB] . the influence 
of the heterogeneity of a coagulated colloidal suspen- 
sion on its mechanical properties was investigated. It 
was shown that more heterogeneous microstructures 
exhibit up to ten times higher elastic properties and 
yield strengths than their more homogeneous counter- 
parts at equal volume fraction. These experiments have 
systematically shown that the local arrangement of the 
powder particles has a strong influence on the mechani- 
cal properties. The relation between the microstructure 
and the mechanical properties, however, is still an open 
question. 

A precondition for any methodical study of the re- 
lation between the microstructure of any granular ma- 
terial and its mechanical properties is the possibility to 
generate particle packings with a controlled degree of 
heterogeneity (DOH). For colloidal particle structures, 
this reproducible control of the heterogeneity is ex- 
perimentally achieved using an enzyme catalyzed gela- 
tion method (DCC = direct coagulation casting [17] 
\TE\). DCC allows for an undisturbed coagulation of 
electrostatically stabilized colloidal suspensions to stiff 
microstructures by an in situ transition of the inter- 
particle potential from repulsive to attractive. This 
can be done along two principal pathways: shifting the 
pH of the suspension to the particles' isoelectric point 
(ApH-mcthod) or increasing the ionic strength in the 
suspension at constant pH (AI- method), which com- 
presses the Debye length of the repulsive inter-particle 
potential. The first method leads to more homogeneous 
microstructures through diffusion-limited aggregation. 
The second method produces more heterogeneous mi- 
crostructures via reaction-rate- limited aggregation [IS] . 

A variation of the ApH-method producing mi- 
crostructures with higher degree of heterogeneity con- 
sists in admixing alkali-swellable polymer (ASP) par- 
ticles to the powder particles under acidic conditions. 
The initial diameter of the ASP particles is small (ap- 
proximately 80 nm) in comparison to the diameter of 
the structural particles (400 nm). The ASP particles 
swell upon increasing pH during the internal gelling 
reaction and unfold to roughly 700 nm in diameter, 
thereby pushing the structural particles in their vicin- 
ity, creating larger pores and thus producing more het- 
erogeneous microstructures. 

One way to investigate the microstructure- 
dependent mechanical properties of coagulated colloidal 
structures is by computational means, such as the dis- 
crete element method or molecular dynamics. These 
methods take intrinsically account of the colloid's par- 
ticulate nature but require to be provided with the ini- 
tial particle configurations. In the case of gravitation- 
ally stable random particle packings with volume frac- 



tions between random loose (<1>rlp ~ 0.55 [2D]) and 
random close packing ('Prop ~ 0.64 [213 ) a large variety 
of algorithms exists. In [22] . for example, an overview of 
available approaches is given. Basically, they can be cat- 
egorized into two main branches: dynamic techniques 
and constructive algorithms. Dynamic techniques, to 
which the method used in this study belongs, rely on a 
gradual swelling of the particles or shrinking of the sim- 
ulation box. Constructive methods, such as the method 
introduced in [22] , construct densely packed particle as- 
semblies using geometrical calculations. Volume frac- 
tions below <Prlp were achieved by molecular dynam- 
ics simulations using long-range attractive interactions 
between the particles and high friction coefficients [23] . 

Despite the large amount of structure generation 
methods, to our knowledge, none of them has been 
shown to provide the possibility to control the mi- 
crostructural arrangement of the particles indepen- 
dently of other parameters such as the volume fraction 
or the particle size distribution. This control, however, 
is crucial to our study of the microstructure-dependent 
mechanical properties and was a motivation to develop 
the void expansion method (VEM) presented in |24) . 
The method was originally inspired by the generation 
of heterogeneous microstructures using ASP particles 
and was shown to allow for a fast and efficient compu- 
tational generation of porous particle structures over a 
broad range of volume fractions. 

In this paper, particle packings with volume frac- 
tions between 0.4 and 0.6 generated using the void ex- 
pansion method are further analyzed in terms of their 
degree of heterogeneity (DOH). The DOH, introduced 
in |25j . represents an abstract concept quantitatively 
describing the heterogeneity of a particle arrangement 
by scalar measures. In order to quantify the DOH, three 
structure characterization methods in combination with 
parameters in fit functions or integral measures are 
used: the pore size distribution, the Voronoi volume dis- 
tribution and the density-fluctuation method. For vol- 
ume fractions between 0.4 and 0.55, we show that the 
number of void-particles strongly influences the DOH 
of the final microstructure. Additionally, for a volume 
fraction of 0.4, the degrees of heterogeneity of struc- 
tures generated using VEM and previous Brownian dy- 
namics (BD) simulations 26,27 a are compared. The 
BD microstructures represent coagulated colloidal mi- 
crostructures, for which the DOH was shown to depend 
on the presence and depth of a secondary minimum in 
the inter-particle potential described by the Derjaguin- 
Landau-Verweg-Overbeek (DLVO) theory [28]. In con- 
trast to these BD-simulations that simulate the coag- 
ulation using widely accepted physical laws and theo- 
ries, VEM is a stochastic method that has the advan- 
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tage of its computational efficiency. We found that the 
DOH-ranges of the VEM- and the BD-microstructures 
overlap to a large extent with the range of the VEM- 
microstructures shifted to slightly higher values. This 
indicates that VEM is well-suited to reproduce coagu- 
lated colloidal microstructures with varying degree of 
heterogeneity. 

2 Materials and Methods 

2.1 Void Expansion Method 

The void expansion method is implemented using the 
particle flow code in three dimensions (PFC 3D ) from 
Itasca Consulting Group, Inc., Minneapolis, Minnesota, 
USA 29 . PFC 3D is based on the discrete element 
method [30) , which allows a modelling of the movement 
of assemblies of rigid spherical particles. In PFC 3D , a 
central-difference scheme is used to numerically inte- 
grate the accelerations and velocities of the particles 
and thus to determine their dynamic behavior. The 
forces on the particles included in our model arise from 
a linear elastic contact law between the particles and 
damping. The contact law is characterized by the par- 
ticles' normal and shear stiffness k n and k s , respectively 
and the damping force is adjusted via the damping co- 
efficient d [31]. 

The void expansion method relies on two kinds of 
particles: "structural particles" and "void-particles". 
The structural particles, with normal stiffness k r ^s and 
shear stiffness fc s ,s, constitute the final microstructures, 
whereas the void-particles (normal and shear stiffness 
k n y and k s y, respectively) are only used during struc- 
ture generation. The edge length Lb ox of the cubic sim- 
ulation box with periodic boundary conditions is calcu- 
lated using the number of structural particles JVg, their 
radius r$ and their volume fraction <&s via 

The Ns structural particles are randomly placed 
in the simulation box. In PFC 3D , particles are not al- 
lowed to overlap during their generation. Thus, in order 
to achieve volume fractions higher than approximately 
0.35, the structural particles are generated with a re- 
duced initial radius of / (m+1), m = 10 in our simula- 
tions. The Ny void-particles having a radius ry *C r$ 
are randomly added to the structural particles. After 
generation of all the particles, the radius of the struc- 
tural particles is increased by means of m repeated ra- 
dius blow-up steps. At each step, the initial particle 
radius is added to the current radius, followed by an 



equilibration of the structure until, after m steps, the 
final particle radius r$ is reached. 

After the structural particles have reached their fi- 
nal size, the radius of the void-particles is cyclically in- 
creased by adding their initial radius ry to their current 
radius in alternation with the performance of 20000 cal- 
culation steps in order to allow for a relaxation of the 
structure. This cyclic increase of the void-particle ra- 
dius simulates the swelling of the ASP particles in the 
experiment and causes the structural particles to rear- 
range and to get in contact with each other. 

The iteration is generally done until the structural 
and void-particles are densely packed and any further 
increase in the void-particle radius leads to a com- 
paction of the structural particles. This is reflected by 
an increasing internal strain energy, which essentially 
depends on the stiffness of the void- and the structural 
particles. In this paper, VEM-microstructures having a 
specific average coordination number are chosen for fur- 
ther analysis. For a volume fraction of 0.4, an average 
coordination number of CV0.4 = 4.7 is targeted. This 
value corresponds to the average coordination num- 
ber of the microstructures resulting from the aforemen- 
tioned BD-simulations, thus allowing for a direct com- 
parison of their degrees of heterogeneity. For volume 
fractions above 0.4, the according coordination num- 
bers are determined using a linear interpolation be- 
tween CV0.4 and CNrcp — 6 [32]. All analyses are 
performed on the structural particles alone, after dele- 
tion of the void-particles. As shown in [23], the evo- 
lution of CN as a function of the void-particle radius 
depends on the void-particle number Ny. The range of 
the ATy -values that allow reaching the targeted CN is 
discussed in Sec. 13.21 in terms of particle overlaps. 

The simulation parameters are summarized in Ta- 
ble [T] In order to reduce the inertia of the void-particles 
their density py was set ten times smaller than the den- 
sity of the structural particles ps, for which the density 
of bulk alumina was chosen. The inter-particle friction 
coefficient p was set to zero in order not to impede any 
particle rearrangements. 



2.2 Quantification of the Degree of Heterogeneity 

Three distinct structural characterization methods are 
used to analyze and assess the heterogeneity of the 
VEM-microstructures [25] : the pore size distribution, 
the density-fluctuation method and the distribution of 
Voronoi volumes. Each method, in conjunction with fit 
functions or integral measures, provides a scalar mea- 
sure that captures and quantifies the DOH of a mi- 
crostructure. 
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Table 1 Simulation parameters 



Parameter 


Symbol 


Value 


Number of particles 


N S 


8000 


Particle radius 


rs 


2.5 x 10" 7 m 


Normal structural particle stiffness 


k n ,s 


10 3 N/m 


Shear structural particle stiffness 


k s ,s 


10" 2 N/m 


Number of void-particles 


N v 


1000 - 16000 


Normal void-particle stiffness 


kn,V 


10 2 N/m 


Shear void-particle stiffness 


k s ,v 


10~ 2 N/m 


Damping coefficient 


d 


0.7 


Friction coefficient 


/' 


0.0 


Volume fraction 


<2> s 


0.4 - 0.6 


Structural particle density 


PS 


3690 kg/m 3 


Void-particle density 


Pv 


369 kg/m 3 



2.2.1 Pore Size Distribution 

The pore size distribution is calculated using the exclu- 
sion probability Ey {r) [33] . It is defined as the probabil- 
ity of inserting a "test" particle of radius r at some arbi- 
trary position in the pore space of a microstructure and 
is estimated using a Monte Carlo approach. The proba- 
bility P{rp > r) of finding a pore with radius rp larger 
than r is obtained using P(rp > r) — J2 r '>r P Ev(r')- 
For coagulated colloidal suspensions resulting from BD- 
simulations it was shown that P(rp > r) follows a com- 
plementary error function given by 

V/r - b 



P{rp > r) — 1 — erf 



ay/2 



(2) 



with a and b the standard deviation and the mean value, 
respectively. Parameter a was shown to nicely reflect 
the structure's degree of heterogeneity, where increasing 
values of a indicate an increasing DOH. 

2.2.2 Density Fluctuation Method 

The density-fluctuation method statistically analyzes 
the spatial distribution of the particle centers as a func- 
tion of grid spacing. Therefor, the cubic simulation box 
is subdivided into n\ cells where n c — 2, . . . ,n™ ax 
{n™ ax = 33, cf. [25] ). The standard deviation a ppc nor- 
malized by the average particle number E ppc as a func- 
tion of grid spacing was shown to depend on the DOH 
of the microstructure and the integral Idf over these 
curves (Eq. ([3])) provides a measure of the DOH. 

tpc (^c) 



'df 



E 



E 



(3) 



ppc 



2.2.3 Voronoi Volume Distribution 



The Voronoi volume of a particle is the volume given by 
the union of all points in pore space closer to the sur- 
face of this particle than to any other particle [33] . The 



Voronoi tessellation thus divides a set of particles into 
a set of space-filling, non-overlapping and convex poly- 
hedrons. In this study, the polyhedron volumes are de- 
termined using the Qhull package [35] . Following a sta- 
tistical mechanics approach, Aste and Di Matteo have 
shown in [36! that the Voronoi volume distribution of 
a particle set follows a so-called fc-gamma distribution 
given by 



f(V f ,k) = 



k k (v f ) k - 
m {vfy 



■ exp(- 



vt 



(4) 



The Voronoi free volume V* = V — V m in is the 
difference between a particle's Voronoi volume V and 
the minimum volume of a Voronoi cell V m i n , which is 
achieved for a regular close packing and is given by 
Vmin = 1-325 V sphere , with V sphe re the volume of a par- 
ticle. The mean Voronoi free volume is a scaling 
parameter and the free parameter k characterizes the 
shape of the curve that very sensitively depends on the 
structural organization of the particles. In (25) . the fc- 
gamma distribution was used to fit the Voronoi vol- 
ume distribution of BD-microstructures. A very good 
agreement between data and fit was obtained and it 
was shown that parameter k reflects very sensitively 
the DOH of the various microstructures. 



3 Results and Discussion 

A qualitative impression of the influence of the void- 
particle number Ny on the resulting microstructures 
is given in Fig. [T] presenting slices with a thickness of 
three particle layers and equal volume fraction of 0.4. 
In the upper row, slices through VEM-microstructures 
generated using N v = 16000 (left) and N v = 1000 
(right), respectively, are shown. In the left slice, the par- 
ticles are rather uniformly distributed, whereas the par- 
ticles in the slice on the right are locally more densely 
packed and thus larger voids can be observed. Very sim- 
ilar properties are observed for the BD-microstructures 
inserted in the lower row. On the left, the most ho- 
mogeneous (surface potential = mV) and, on 
the right, the most heterogeneous BD-microstructure 
(<Pb = 15 mV) are shown. Figure Q] suggests that the 
use of larger void-particle numbers leads to more ho- 
mogeneous microstructures and that VEM allows ob- 
taining very similar microstructures as BD in terms of 
heterogeneity, as supported by the quantitative analy- 
ses in Sec. 13.11 In Sec. 13.21 the DOH of the microstruc- 
tures generated using VEM and BD are compared and, 
in Sec. 13.31 the possibility to generate microstructures 
with a varying DOH for volume fraction above 0.4 is 
investigated. 
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Fig. 1 Slices through VEM-microstructures (upper row) with a 
more homogeneous microstructure (Ny = 16000, left) and a more 
heterogeneous microstructure (Ny = 1000, right). The lower row 
shows the most homogeneous ('Z'o = mV, left) and most het- 
erogeneous microstructure (<?o = 15 mV, right) resulting from 
earlier BD-simulations 26 . Slice thickness: three particle diam- 
eters; particle diameter = 0.5 /xm, volume fraction = 0.4 (for all 
structures). 

3.1 Influence of the Void-Particle Number 

In this section, the results obtained for the three dis- 
tinct methods used to quantify the DOH are presented 
in detail for the VEM-microstructures with a volume 
fraction of 0.4 and varying void-particle number Ny. 

3.1.1 Pore Size Distribution 

The probability P(rp > r) of finding pores with a ra- 
dius rp larger than r is shown in Fig. [5] for various 
VEM-microstructures with void-particle numbers rang- 
ing from 1000 to 16000 (symbols). The data is shown 
as a function of r normalized by the particle radius 
rs- For a given pore radius rp > 0, P(rp > r) de- 
creases for increasing void-particle numbers Ny. The 
probability of finding pores with a radius larger than 
0.5 rs is 1.6 times higher in the Ny = 1000 than in the 
Ny — 16000 microstructure. Finding pores larger than 
0.75 r 5 and 1.0 rs is 4.3 and 18.7 times more probable 
in the microstructure with Ny = 1000 than in the one 
with Ny = 16000, respectively. The solid lines in Fig. [2] 
represent the complementary error function fit curves 
obtained using Eq. @. Fit parameters a and b and the 
corresponding R 2 ^-values are summarized in Table [3J 
The R a b - values close to one indicate excellent fits. The 



Table 2 Measures of the degree of heterogeneity as a function of 
the void-particle number Ny. fit parameters a and k, correspond- 
ing to the widths of the pore size and Voronoi volume distribu- 
tion, respectively, and the integral over the density fluctuation 
curves Im. Additionally, the second fit parameter of the cumu- 
lative pore size distribution (b) and the corresponding i? 2 -values 
are shown. 



N v 


o 


b (io- 2 ) 


a,b 




k 


K 


1000 


0.605 


-7.44 


0.997 


24.2 


2.46 


0.977 


2000 


0.530 


-4.63 


0.999 


23.6 


2.99 


0.995 


4000 


0.472 


-2.56 


0.999 


23.1 


3.61 


0.998 


8000 


0.431 


-1.21 


0.999 


22.7 


4.67 


0.999 


11000 


0.416 


-0.652 


0.999 


22.6 


5.02 


0.999 


13000 


0.411 


-0.458 


0.998 


22.6 


5.50 


0.999 


16000 


0.401 


-0.0227 


0.999 


22.5 


5.74 


0.998 
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Fig. 2 Probability P(rp > r) of finding pores with a radius rp 
larger than r vs. r/rg obtained using the pore size distribution 
for the various VEM-microstructures with <Pg = 0.4 (symbols). 
Solid lines denote the corresponding fits using a complementary 
error function (Eg. [2t. 

largest value a is found for the microstructure gener- 
ated using 1000 void-particles. For increasing values of 
Ny, parameter a decreases, which reflects a decreasing 
probability of finding larger pores and thus a decreasing 
DOH. 



3.1.2 Density Fluctuation Method 

The density fluctuations of the various VEM- 
microstructures are presented in Fig. [3] as a function 
of the grid spacing normalized by the particle diameter 
ds = 2r$. The density fluctuations decrease with in- 
creasing void-particle number. This behavior is quanti- 
fied using Idf as given in Eq. ([3]). The values are summa- 
rized in Table [21 Idf increases with decreasing Ny and 
thereby reflects an increasingly heterogeneous reparti- 
tion of the particles. 
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Fig. 3 Density fluctuations as a function of grid spacing for the 
various VEM-microstructurcs with <P$ = 0.4. 
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Fig. 4 Voronoi volume distribution P(u) (symbols) as a function 
of v and corresponding t-gamma fits (lines) for various VEM- 
microstructures (#g = 0.4). 



3.1.3 Voronoi Volume Distribution 

The Voronoi volume distribution P(y) of the vari- 
ous VEM-microstructures is shown in Fig. [4] for void- 
particle numbers ranging between 1000 and 16000 
(symbols) as a function of v = —j, the free Voronoi 
volume normalized by the mean free volume. For in- 
creasing Ny, the peak height increases and, conse- 
quently, the width of the curve decreases. This behavior 
indicates that more homogeneous microstructures are 
found toward increasing Ny, which is confirmed by a 
fit of the curves using the /c-gamma distribution given 
in Eq. Q. The fits are shown as lines in Fig. [4] The cor- 
responding values for the parameter k, summarized in 
Table [H increase for increasing Ny reflecting the shift 
toward more homogeneous microstructures. Very good 
fit results are achieved as indicated by the i?^-values. 
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Fig. 5 Interdependence of the DOH-measures (<Pg = 0.4): In- 
tegral over the density fluctuation curves I^f (circles, left scale) 
and width of the Voronoi volume distribution 1/k (triangles, right 
scale) as a function of the width of the pore size distribution a 
for the various VEM-microstructures (full symbols) and for the 
BD-microstructures (open symbols |25p at = 0.4. 



3.2 Degree of Heterogeneity for VEM- and 
BD-Structures at volume fraction 0.4 

The interdependence between the three measures of the 
degree of heterogeneity is shown in Fig. [5] On the left 
scale, the integral over the density fluctuation curves 
Idf is shown as a function of the width of the pore size 
distribution, given by parameter a. The corresponding 
values are shown as circles. The right scale presents 1/k 
reflecting the width of the Voronoi volume distribution 
in dependence of a, shown as triangles. Full symbols de- 
note the values obtained for the VEM-microstructures. 

The linear arrangement of the two curves sug- 
gests a pairwise affine relation between the various 
DOH. Furthermore, the values found for the VEM- 
microstructures are very close to those obtained for the 
BD-microstructures, which are inserted as open sym- 
bols. These two points substantiate the equivalence of 
the three methods as already observed in [2"5] . 

The DOH-range covered by the VEM- 
microstructures, in terms of parameter a, extends 
from 0.4 to 0.6. This overlaps to a large extent with 
the range of the BD-microstructures from 0.37 to 
0.54. Thus, the range covered by VEM is roughly 
20% broader starting at a slightly higher DOH than 
the range of the BD-microstructures considered in this 
work. 

The sensitivity of the DOH on Ny is demonstrated 
in Fig. [6] presenting 1/k relative to its maximum value 
(l/k) max , achieved for N v = 1000, as a function ofN v . 
After a strong initial decrease for small values of Ny, 
1/k levels off toward higher values of Ny- A further 
increase of Ny has only a negligible influence on the 
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q ^ most homogeneous BD-microstructure (4^ = mV)_ 

4000 8000 12000 16000 20000 

N v 

Fig. 6 DOH-mcasure 1/fe relative to (l/k) max as a function 
of the void-particle number Ny for <Ps = 0-4 (symbols). The 
dashed lines denote the corresponding values of the most and 
least heterogeneous BD-microstructure 25 . 

DOH as shown by the insertion of the additional 1/fe- 
value for Ny = 20000. Indeed, an increase in Ny from 
16000 to 20000 decreases the DOH by only 2%, which 
suggests that the DOH of the most homogeneous BD- 
microstructure, represented by the lower dashed line in 
Fig. El cannot be reached using VEM. This is supported 
by the fit of the data using 

(1/fe - 1/fco) oc exp (j3N v ) , (5) 

yielding a minimum DOH of 1/feo = 0.43, which is 
larger than the 1/fe- value of the homogeneous BD- 
microstructure. The second fit parameter and the cor- 
responding correlation coefficient are (3 — —0.25 x 10~ 3 
and R 2 = 0.99, respectively. The reason for this asymp- 
totic behavior is most probably related to the VEM- 
algorithm itself in combination with the chosen targeted 
mean coordination number CN = 4.7. As shown in [24] , 
CN as a function of the void-particle radius follows a 
characteristic step-like shape with CNi the coordina- 
tion number at the curve's inflection point. CN was 
found to decrease for increasing void-particle numbers 
Ny- Above the inflection point, CN scales as a power 
law with exponent 0.37, independently of Ny- There- 
fore it becomes increasingly difficult to reach a tar- 
geted CN > CN toward larger values of Ny, at least 
without accepting considerable particle overlaps and 
consequently high internal stresses. The mean particle 
overlap, as determined by the maximum in the pair- 
correlation function was found to be 1.4% of the particle 
diameter d s in the case of Ny = 13000 and CN = 4.7. 
For Ny = 16000 the mean overlap is 1.5% ds, how- 
ever with a slightly lower CN of 4.63. The maximum 
particle overlap is below 2% ds for all microstructure 
with Ny < 8000 and increases to roughly 4% ds for 
Ny — 16000. For Ny = 20000 a mean and maximum 



overlap of 2% ds and 5% ds, respectively was already 
found for CN = 4.57, which led to the conclusion that 
Ny — 16000 and thus a void- to structural particle 
number ration n max = 2.0 constitutes the upper limit 
of VEM, for a volume fraction of 4>s = 0.4. 

The lower limit of Ny, at around 1000, is mainly 
imposed by the decreasing fit quality obtained in the 
case of the pore size and Voronoi volume distribution, 
as expressed in terms of the i? 2 -values (Table [2]). The 
minimum void- to structural particle number ratio in 
the case of <P S = 0.4 is thus n min = 0.125. 

3.3 Influence of the Volume Fraction on the DOH 

In the following, the DOH-range is investigated for 
volume fractions above 0.4. In particular, VEM- 
microstructures with volume fractions <Ps — 0.45, 0.50, 
0.55 and 0.60 have been generated and analyzed anal- 
ogously to the structures with <Ps = 0.4. This range of 
<?5 was chosen since it is of particular interest for fur- 
ther simulations of the microstructure-dependent prop- 
erties of coagulated colloidal suspensions [16 . As men- 
tioned in Sec. 12.11 the targeted average coordination 
numbers for the structures with ^5 > 0.4 are obtained 
using a linear interpolation between CNqa = 4.7 and 
CNrcp = 6.0 yielding approximately 5.0, 5.25, 5.5 and 
5.75, respectively. 

In this section, the DOH is measured by means of 
parameter a obtained by fitting the cumulative pore 
size distribution using the complementary error func- 
tion given in Eq. [2 This method was chosen because it 
yielded the best fit results over the whole range of vol- 
ume fractions above 0.4 (R 2 > 0.999 for all fits). The re- 
sults are compiled in Table[3J For each volume fraction, 
an increasing void-particle number entails a decreas- 
ing parameter a and therefore a decreasing DOH. This 
is summarized in Fig. presenting the DOH-values of 
these structures, as measured by parameter a (circles). 
Particular interest is given to the maximum and min- 
imum values of a as a function of <Ps'- a max (<Ps) and 
a mm (@s), respectively. Our results suggest linear de- 
pendencies as illustrated by the solid lines in Fig. [7] 
The respective fits yield: 

a" UB (#s) = -2.16 $s+ 1-47 

a min ($ s ) = -1.11 $ s + 0.842. [ } 

which, for the DOH-range measured by parameter a, 
results in 

Aa{$ s ) = -1-05 <P S + 0.624. (7) 

Equation [7] describes the narrowing of the DOH-range 
toward increasing volume fraction. With respect to the 
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Table 3 DOH-measure a as a function of the void-particle num- 
ber Ny for various volume fractions <P g > 0.4. R 2 > 0.999 for all 
fits. 
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Fig. 7 DOH-range for volume fractions <Ps between 0.4 and 0.6 
as expressed by measure a obtained from the pore size distribu- 
tion. 

DOH-range at <P S = 0.4, the ranges for ^ s = 0.45, 0.5 
and 0.55 are reduced by 30, 55 and 76%, respectively. 
The linear extrapolations of a max (<P s ) and a mm (<?s) 
toward higher volume fractions are shown as dashed 
lines in Fig. [7j The lines intersect at a volume fraction 
of 0.594 with corresponding DOH- value ao — 0.183. 

The DOH-range for a random packing of monosized 
spheres is thus reduced to a unique value at this vol- 
ume fraction. This was confirmed using two additional 
VEM-microstructure at <Pg = 0.6 and using 500 and 
8000 void particles (triangles in Fig. [7J. The DOH- 
measure a for these two structures yields 0.187 and 
0.191, respectively. The range is thus narrowed down to 
roughly 2% with respect to the DOH-range at <!>$ — 0.4. 

The decreasing DOH-range toward larger volume 
fractions is qualitatively explained by the decrease of 
the pore space volume required for particle rearrange- 
ments. The volume fraction of 0.6, however, is below the 
RCP-limit, which was assumed to be the volume frac- 
tion, at which variations in DOH become zero. In [57] , 
for example, the nature of the RCP limit was studied 
by means of a statistical mechanics approach and it was 



found that the entropy and thus the number of accessi- 
ble configurations of a disordered packing of monosized 
spheres reaches a minimum at the RCP limit. Addition- 
ally, a sharp decrease in parameter k was found as the 
volume fraction crosses the RLP limit in [38 . Following 
our interpretation of k, this sharp decrease corresponds 
to a sharp increase in DOH, which is not found in our 
data. More detailed studies are required in order to in- 
vestigate these differences. 

4 Summary and Conclusions 

In this paper, the degree of heterogeneity of microstruc- 
tures generated using the void expansion method has 
been analyzed and quantified using three distinct tech- 
niques based on the pore size distribution, the density- 
fluctuation method and the Voronoi volume distribu- 
tion. In particular, the influence of the void-particle 
number and the volume fraction on the DOH was in- 
vestigated. 

The various pore size distributions were fitted using 
a complementary error function achieving very good fit 
results over the whole range of volume fractions investi- 
gated in this study. For all volume fractions, parameter 
a, corresponding to the width of the pore size distribu- 
tion, increases for decreasing Ny, thereby reflecting the 
transition toward a higher DOH. For this range of vol- 
ume fractions, the void expansion method thus allows 
for a generation of particle arrangements, for which the 
DOH can be controlled quasi-continuously over a broad 
range by means of the void-particle number. In partic- 
ular, the DOH-range becomes narrower and shifts to 
lower values as a function of increasing volume fraction. 
Our results indicate that the DOH-range is reduced to a 
single value at a volume fraction of roughly 0.6, which, 
interestingly, is below the RCP limit. Further studies 
are required in order to analyze the behavior of the 
DOH and its range for volume fractions toward and 
beyond the RCP limit. 

For VEM-structures with volume fraction 0.4, the 
significant influence of the void-particle number on the 
DOH as measured by parameter a is confirmed by the 
density-fluctuation method and the Voronoi volume dis- 
tribution. I df , given by the integral over the various 
density fluctuation curves, increases for decreasing Ny. 
Parameter k, characterizing the shape of the Voronoi 
volume distributions, decreases for decreasing Ny. 

The degree of heterogeneity of the VEM- 
microstructures with a volume fraction of 0.4 has 
been compared to that of BD-microstructures [26l 
127] analyzed in [35]. The structures generated using 
BD-simulations represent coagulated colloidal particle 
structures where the coagulation was simulated using 
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widely accepted physical laws and theories. The degree 
of heterogeneity of the BD-microstructures was shown 
to be closely related to the presence and depth of a 
secondary minimum in the DLVO potential, which 
is determined by the particles' surface potential. In 
contrast to BD, VEM is a stochastic method that 
mimics the swelling of ASP particles in experiment, 
having the advantage of its computational effective- 
ness. The DOH-range covered by VEM in terms of 
parameter a, characterizing the width of the pore 
size distribution, is approximately 20% broader and 
shifted to slightly higher values than the DOH-range 
of these BD-microstructures. The most homogeneous 
BD-microstructure, for which a surface potential 
•^0 = mV was used, could not be reproduced by 
VEM in terms of heterogeneity. A further decrease 
of the DOH through an increase in Ny is mostly 
impeded by the choice of the targeted coordination 
number (roughly 4.7). For void-particle numbers 
exceeding 16000, this value cannot be attained without 
considerable particle overlap and therefore internal 
stress in the microstructure. On the other hand, VEM 
allows generating more heterogeneous microstructures 
than these BD-simulations, where the increasing depth 
of the secondary minimum toward higher surface 
potentials inhibits a complete coagulation. 

This study has undoubtedly shown, that the void- 
particle number constitutes an important simulation 
parameter of VEM with regard to the DOH of the fi- 
nal microstructure. This parameter indeed allows for a 
quasi-continuous control of the DOH over a broad range 
of volume fractions. In order to compare the structural 
arrangement and to assess the influence of the void- 
particle number on the final microstructures, the con- 
cept of the DOH itself has proven particularly useful. 
The comparison of VEM- and BD-generated structures 
at a volume fraction of 0.4 allows concluding that VEM 
nicely reproduces the microstructures of coagulated col- 
loids. Thereby, the void expansion method facilitates 
the further study of the microstructure-dependent me- 
chanical properties of coagulate colloidal structures or 
granular matter in general. 
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